Dynamical mean-field theory of photoemission spectra of actinide compounds 
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A model of photoemission spectra of actinide compounds is presented. The complete multiplet 
spectrum of a single ion is calculated by exact diagonalization of the two-body Hamiltonian of the 
/" shell. A coupling to auxiliary fermion states models the interaction with a conduction sea. The 
ensuing self-energy function is combined with a band Hamiltonian of the compound, calculated in 
the local-density approximation, to produce a solid state Green's function. The theory is applied to 
PuSe and elemental Am. For PuSe a sharp resonance at the Fermi level arises from mixed valent 
behavior, while several features at larger binding energies can be identified with quantum numbers 
of the atomic system. For Am the ground state is dominated by the |/ 6 ; J = 0) singlet but the 
strong coupling to the conduction electrons mixes in a significant amount of f 7 character. 

PACS numbers: 71.10.-w, 71.27. +a, 71.28.+d, 78.20.Bh 



The electronic structure of actinides has been 
subject of extensive experimentali*2iMi5i& and 
theoretical£i2i±fliiii±2ii3' 14 - 15 investigations in recent 
years. The key issue is the character of the /-electrons 
which show varying degrees of band- like and/or localized 
behavior. Photoemission in particular provides energy 
resolved information about the /-electrons, the under- 
standing of which is far from completed The actinide 
antimonides and chalcogenidesi& display narrow band- 
like features around the Fermi level as well as atomic-like 
features at higher binding energies. Similarly, among 
the elements a distinct shift of /-character away from 
the Fermi level happens from Pu to Amj2*il where also 
the cohesive properties suggest a shift from itinerant to 
localized behavior of the /-electrons. The atomic-like 
features of photoemission spectra are not well analyzed 
in terms of multiplets of isolated actinide atoms, in 
contrast to the situation for rare earth compounds^ On 
the other hand, band structure calculations generally 
lead to narrow /-bands pinned at the Fermi level but 
no structures at large binding energy, and additional 
modelling has been introduced to account for the dual 
character of the /-electrons £2ii2iiiii£ii£i The recent 
development of dynamical mean-field therory (DMFT) 
has in particular spawned advancements 

In the present work a novel model of actinide pho- 
toemission is presented, which is capable of describing 
both high and low energy excitations. In applications 
to PuSe and Am metal it is shown that interaction of 
the /-electrons with the conduction electrons leads to 
complex ground state configurations for the solid, which 
account for the non-trivial features of the photoemission 
spectra. The theory is based on the Hubbard-I approach 
suggested in Ref . [l!| but augmented with coupling to the 
conduction sea. 

To model the photoemission of actinide compounds 
the k-integrated spectral function A(ui) — 7r~ 1 ImG' oc (a;), 
with 



G 1oc {uj) 



-E 



is calculated. A dynamical mean-field, is added 

to a band structure Hamiltonian, which is calculated by 
standard density-functional-theory based methods^ and 
the crystal Green function Gk obtained by inversion: 



Gk(w) 
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is calculated from an effective impurity model de- 
scribing an isolated atom coupled to a bath of uncorre- 
cted electrons: 



Hi, 



H r , 



(3) 



The atomic Hamiltonian, H at0 m, includes the electron- 
electron interaction, V ee , within the /" shell of an isolated 
actinide atom: 



H„ 



(1) 
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Here, indices rrii = 1, .., 14 refer to the individual orbitals 
of the /" shell including spin, and /+ and f m are the cre- 
ation and annihililation operators. The bare /-level is ej, 
and [i is the chemical potential. The second term is the 
spin-orbit energy, £ being the spin-orbit constant, which 
is calculated from the self-consistent band structure po- 
tential, and Si and h are the spin and orbital moment 
operators for the t'th electron. 

The two-electron integrals of the Coulomb operator, 
U mim2Tn3mi = (mim2\V ee \m 3 m4} , may be expressed in 
terms of Slater integrals, F k , k — 0,2,4,6^ which 
are computed with the self-consistent radial /-waves 
of the band structure calculation. In practice we find 
it necessary to reduce the leading Coulomb integral 
F° = U„ lim2mi m 2 from its bare value, similar to other 
studies»2i The physical reason for this is the significant 
screening of the charge-fluctuations on the /-shell by the 
fast conduction electrons, which happen as part of the 
photoemission process, but which is not treated in the 
present model. A complete treatment would involve the 
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computation of the solid dielectric function, as is done in 
the GW-approximationi^2i£i The higher Slater integrals, 
F 2 , F A and F 6 , govern the energetics of orbital fluctua- 
tions within the f n shell, for which the screening is less 
important, and we therefore use the ab-initio calculated 
values for those parameters. 

The coupling to a sea of conduction electrons is mod- 
elled as a simple hopping term for each of the 14 /- 
electrons into an auxiliary state at the Fermi level: 



(5) 



The impurity Hamiltonian, Eq. Q , is solved by exact di- 
agonalization, and the dynamical mean-field, S mm /, fol- 
lows from the impurity Green's function, (G lmp ) mm i , 
which is computed from the eigenvalues, E^, and eigen- 

states, of Himp as: 

lu + &u — hj v 

fJ,V p 



In the zero-temperature limit, \ 
or \v) is the ground state, 
wise. Further, G° = (uj — e/ 
function of the bare /-level. 



= 1 



if either state 
and = other- 

+ it) 1 is the Green's 
Alternatively, G° may 



be determined by the DMFT self-consistency through 
(G )- 1 = (G 100 )" 1 + E(lj), whith G loc given by Eq. (1). 
In the present applications there are only insignificant 
difference between these two expressions for G°, i.e. the 
DMFT self-consistent solid state local Green's function 
is fairly close to that of a bare f-electron (at the energy 
position given by the LDA eigenvalue) . The effect of the 
coupling term in Eq. (|SJ| is to mix the multiplet eigen- 
states corresponding to fluctuations in /-occupancy. The 
eigenstates of Hi mp will not contain an integral number of 
/-electrons, but may deviate more or less from the ideal 
atomic limit, depending on parameters, in particular the 
coupling strength V. 

In a more realistic treatment the conduction states in 
(JSJ) should aquire a width, and the hybridization parame- 
ter V be energy dependent^ The energy dependent hy- 
bridization may be calculated directly from a LDA band 
structure calculation (or from the self-consistent DMFT 
bath Green's function)^ but the straightforward map- 
ping to the above model leads to overestimation of the in- 
teraction, since the auxiliary electrons cost no excitation 
energy. Hence, V is treated as an adjustable parameter 
of the present theory. The hopping parameter V in (JSJ) 
could in principle depend on orbital index (which would 
be relevant in highly anisotropic crystal structures^), 
but here we apply the isotropic model of Eq. (JSJ . Simi- 
larly, with no particular extra effort the bare /-energies in 
(@J could be orbital dependent, e.g. for studies of crystal 
field effects. In the present study only a single /-energy 
is used, given by its average value in the band structure 
calculation. 
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FIG. 1: Illustration of the lowest multi plet s for a Pu atom 
at parameter values pertinent to PuSe |24|. The energy is 
reckoned from the lowest |/ 5 ; J = 5/2} state. Only the lowest 
few f n states are shown. Each level is labeled by its (approx- 
imate) Russell- Saunders term. The ground state, when the 
coupling to the conduction states is switched on, is formed 
primarily as a linear combination of the lowest lying / 5 and 
/ 6 states. The energy gain due to the coupling is 0.14 eV; for 
clarity the effect is exaggerated on the figure. 



The theory outlined above may be applied to any /- 
element, i.e. any /-occupancy. The chemical potential 
fj, determines the ground state of the model. Thus, this 
parameter must be adjusted to match the /-element un- 
der study. The adjustment of /z to a certain extent bal- 
ances any inaccuracies in the values of the primary Slater 
integral, F°, since a shift in the latter may be compen- 
sated by a shift in fi. With F° and /i fixed, the ener- 
gies of the ground and excited state multiplets of H atom , 
\f n ; JM) are determined. Fig. 1 ilustrates this for pa- 
rameters relevant for a Pu atom in PuSeiSiS^ This se- 
lects |/ 5 ; J = 5/2) as the ground state with |/ 6 ; J = 0) 
at marginally 44 meV higher energy. 

Introducing the coupling to the auxiliary states, Eq. 
©, leads to a mixed ground state, composed of 70 % 
f and 30 % / 6 , with 67 % and 24 % stemming from 
the lowest |/ 5 ;J = 5/2) and |/ 6 ; J — 0) states, respec- 
tively. Similar mixing occurs in all excited states, and 
all the degeneracies of the atomic multiplets are com- 
pletely lifted by the coupling. In particular, the ground 
state is non-magnetic, and the local-moment character 
of the / 5 ion will only commence at elevated tempera- 
tures. A number of eigenstates occur at low excitation 
energies and allow for low energy transitions in the pho- 
toemission. These transitions would be absent in a pure 
|/ 5 ; J = 5/2) ground state. 

The calculated spectral function, Eq. Q), of PuSe is 
depicted in Figure 2 and compared to the experimental 
photoemission spectrum of Ref. Q. There is a one-to- 
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FIG. 2: (Color online) Calculated spectral function for PuSe, 
/-part in full line (red) and non-/ in dashed line, compared 
to the experimental photoemission spectrum of Ref. [I] (in 
dash-dotted and blue) at 40.8 eV photon energy. The energy 
is measured relative to the Fermi level, 1261) - The dominating 
final state atomic term is indicated above each major peak. 
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FIG. 3: (Color online) Calculated spectral function for ameri- 
cium metal, /-part in full line (red) and non-/ in dashed line, 
compared to the experimental photoemission spectrum of Ref. 
ITU (in dash-dotted and blue) recorded at 40.8 eV photon en- 
ergy. The energy is measured relative to the Fermi level, |2fi| . 
The dominating final state atomic term is indicated. 



one correspondence between the measured and calculated 
spectral features, and the theory allows identification of 
each of these. In particular, a sharp resonance at the 
Fermi level is due to the low-energy 7 Fo — > 6 H 5 / 2 excita- 
tions, which may be identified with heavy-fermion type 
behavior £ At 0.4 eV binding energy excitations into the 
6 H 7 / 2 final state give rise to a small peak, while the more 
pronounced peak at ~ 1 eV binding energy marks the 
excitations into the 6 F 5 / 2 final state, which also has a 
spin-orbit satellite, visible around 1.4 eV binding energy. 
The major emission peak around 2 eV binding energies 
is due to 6 H 5 / 2 — > 5 l4 transitions with small 7 Fo — >6 Ps/2 
contribution. At even larger binding energy a broad emis- 
sion coincides with the Se p-band, (dashed line in Fig. 2). 
It is unclear whether the broad shoulder seen in the ex- 
periment between 4 and 5 eV includes emission related 
to /-features or is solely due to the Se p-band. The iden- 
tification of PuSe as a mixed- valent compound was made 
by Wachter^ and the present work supports this view. 
Details of the identifications of the spectral features are 
different, though, in particular in the present theory it 
is not necessary to reduce the multiplet splittings within 
each /" configuration compared to the atomic values. 
The matrix elements in the numerator of Eq. (6) en- 
sures that all atomic selection rules are obeyed, and de- 
termine the weights of the individual peaks. Left out of 
the present treatment is the matrix element between the 
/-electron one-particle wave and the the photoelectron, 
which may be assumed to vary only slowly with energy 
over the range of the valence bands. However, the exper- 
imental spectrum contains contributions from all allowed 
transitions, /- as well as non- /-related. The justification 
for the comparison in Fig. 2 (and in Fig. 3 for ameri- 
cium below) of the theoretical /-spectral function with 
the experimental photoemission spectrum recorded with 
40.8 eV photons is that the cross section strongly favors 
emission of an /-electron at this energy. The experiments 
reported in Ref. Q include photoemission spectra taken 
with 21.2 eV electrons, at which energy the contribution 
of the non-/ electrons dominates, with a spectrum quite 
close to the non-/ spectral function shown in Fig. 2. Fi- 
nally, life-time effects and secondary electron processes 
are not considered in the present theory. 

The calculated spectral functional of Am is compared 
to photoemission experiment' 17 in figure 3. The multi- 
plet level diagram of the isolated Am atom differs from 
that of Pu shown in Fig. 1 by shifting the balance 
such that |/ 6 ; J = 0) is now the ground state with the 
|/ 7 ; 8 S 7/2 ) at 0.7 eV higher energy, and \f 5 ; 6 H 5/2 ) at 
2.3 eV higher energy. Spin-orbit interaction seriously vi- 
olates the Russell-Saunders coupling scheme. The overlap 
of the americium ground state and the Russell-Saunders 
|/ 6 ; 7 Fq) state is 0.69. The influence on the photoemis- 
sion transition probabilities is substantial. For example, 
the low-energy transitions from the |/ 6 ; J — 0) ground 
state to the \f 5 ; 6 H 5 / 2 ) and \f 5 ', 6 H 7 / 2 ) occur in ratio 
2.6 : 1, while Russell-Saunders coupling predicts the ratio 
1 : 2.5ji& i.e. exactly reversed weight on the two compo- 
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FIG. 4: (Color online) Calculated /-part of the spectral func- 
tion for americium metal, for varying strength of the coupling 
to the conduction electrons, a) V = 0.0 eV , and b) V = 0.16 
eV (dash-dotted and blue), V — 0.33 eV (full line and red), 
and V — 0.49 eV (dashed and green). The spectrum pre- 
sented in Figure 3 corresponds to the full line curve of the 
present plot with a larger broadening |26|. 



nents. For the same reason, the J = 7/2 satellite of the 
peak at the Fermi level in PuSe, Fig. 1, is significantly 
weaker compared to the main peak. 

The coupling to the conduction electrons, Eq. JSJ, 
leads to further discrepancy with the free atom Russell- 
Saunders ground state. Good agreement with the experi- 
mental spectrumii is obtained with a coupling parameter 
V = 0.33 eV in Eq. JSJ. The ensuing Am ground state 
remains a singlet but attains a significant (16 %) admix- 
ture of |/ 7 ; 8 SV/2) character. The |/ 6 ; J = 0) state has 
still got 61 % weight in the coupled ground state, while 



the remaining weight scatters over the higher levels of 
the / 6 configuration. The effect on the spectral function 
is drastic, in particular by leading to significant emis- 
sion of |/ 7 ; 8 <SV/ 2 ) — ► \f e '-J) peaking at 1.7 eV binding 
energy with a broad tail towards lower binding energy, 
i.e. this transmission accounts for the shoulder observed 
in the experiment around 1.8 eV. The main emission of 
|/ 6 ; 7 Fo) -> If-^H^F^P) occurs from - 2.5 eV and 
higher binding energies, and accounts for the plateau seen 
in both experiment and theory between 2.5 and 3.5 eV. 

In Figure 4 the evolution of the calculated Am spec- 
trum with the strength of the coupling to the conduction 
electrons is shown. In Fig. 4(a) the coupling has been set 
to zero, and the spectrum shows the dominating emission 
between 2 and 4 eV binding energy due to |/ 5 ; 6 #5/2) and 
|/ 5 ; 6 F 5 / 2 ) final states and their J — 7/2 satellites. The 
small peaks at -5.1 and -6.2 eV are remnants of the emis- 
sion to |/ 5 ; 6 -P5/2) final states, highly distorted due to 
spin-orbit interaction in both initial and final state. In 
Fig. 4(b) the distortion of the spectral function due to 
coupling to the conduction sea is illustrated, at V = h, 
1 and | times the value used in Figure 3. At binding 
energies less than 1.7 eV the triangular structure evolves 
due to |/ 7 ; 8 67/2} admixture into the initial state. At the 
highest value of V the emission broadens considerably. 
At the Fermi level a peak grows up with increasing cou- 
pling strength, visible as a small shoulder in Figure 3, 
but not resolved in the experimental data. 

In conclusion, a theory has been developed which al- 
lows a detailed description of the photoemission spectra 
of actinide compounds. The three most important en- 
ergy scales are the / intra-shell Coulomb interactions, the 
spin-orbit interaction, and the coupling to the conduction 
electrons, which are all incorporated in the theory. All 
atomic selection rules for the photoexcitation are obeyed. 
The key approximation is that of a dynamical mean-field 
calculated for a single actinide ion with a simplified in- 
teraction with conduction electrons. Compared to recent 
advancements in the field 8 ^^ the present approach treats 
the impurity atom essentially exact including interaction 
with the conduction sea. Hence, the theory applies best 
to systems on the localized side of the Mott transition of 
the / shell. In particular, for PuSe and Am metal the 
interaction leads to formation of complex ground states 
which strongly influence the photoemission spectra, since 
the initial states of the photoemission process are not 
free-atom like. The embedment into the solid via the 
LDA band Hamiltonian in Eq. maintains and broad- 
ens the atomic features. The DMFT self-consistency cy- 
cle has only minute influence on the calculated spectra. 
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electron'. Support from the Danish Center for Scientific 
Computing is acknowledged. 
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